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ABSTRACT 

Context. The weak orbital-phase dependent reflection signal of an exoplanet contains information on the planet surface, such as the 
distribution of continents and oceans on terrestrial planets. This light curve is usually studied in the time domain, but because the 
signal from a stationary surface is (quasi)periodic, analysis of the Fourier series may provide an alternative, complementary approach. 
Aims. We study Fourier spectra from reflected light curves for geometrically simple configurations. Depending on its atmospheric 
properties, a rotating planet in the habitable zone could have circular polar ice caps. Tidally locked planets, on the other hand, may 
have symmetric circular oceans facing the star. These cases are interesting because the high-albedo contrast at the sharp edges of the 
ice-sheets and the glint from the host star in the ocean may produce recognizable light curves with orbital periodicity, which could 
also be interpreted in the Fourier domain. 

Methods. We derive a simple general expression for the Fourier coefficients of a quasiperiodic light curve in terms of the albedo 
map of a Lambertian planet surface. Analytic expressions for light curves and their spectra are calculated for idealized situations, and 
dependence of spectral peaks on the key parameters inclination, obliquity, and cap size is studied. 

Results. The ice-scattering and ocean glint contributions can be separated out, because the coefficients for glint are all positive, 
whereas ice sheets lead to even-numbered, higher harmonics. An in-view polar cap on a planet without axial tilt only produces a 
single peak. The special situation of edge-on observation, which is important for planets in transit, leads to the most pronounced 
spectral behavior. Then the respective spectra from planets with a circumventing ocean, a circular ocean (eyeball world), polar caps, 
and rings, have characteristic power-law tails rT^', and (-l)"^*n“^. 

Conclusions. Promising recently discovered planetary systems may be selected as candidates for long-term (multiyear) observation: 
their Fourier spectra could separate the different planets and reveal or identify a water-covered planet with polar caps. 

Key words. Planets and satellites: detection - surfaces - oceans - rings - Methods: analytical -Techniques: photometric 


1. Introduction 


In the past decade more than a thousand planets near other stars 
have been found, proving that nearly every star has planets. The 
methods of discovery range from direct detection, by watching 


their birth in the circumstellar planetary disks (Greaves et al. 
|2005|l, to more indirect detection in the stellar light dip during 


the passage of a transiting planet in front of its host star (|Hidas 


jet al.||20()5] Steffen et al.|2012 1 and in the time-varying phases 


from reflected star light |Gaulme et al.||2010[ ). Indirectly infer¬ 
ring the planet’s existence is possible by monitoring the motion 
of the host star ( [Benedict et al.|20M] l about the common center 
of mass via measurement of the star’s position coordinates or its 
Doppler-shifted velocity ( Mayor & Queloz|1995 1 , or through the 
detection of stellar light variations due to the effect of the ellip¬ 
soidal deformation of the star by a planet ([Welsh et al.|2010|). 


The phase light curve, which is the (for a twin Earth 10 
times) weaker overall reflection signal of an orbiting planet on 
top of the direct stellar light, is a periodic function in time. This 
periodicity arises as the following combination of different ef¬ 
fects: (i) Eor a noncircular orbit the distance to the star varies, 
so the planet generally scatters more light at periastron and less 
at a-periastron (see Appendix), (ii) Eor inclined observation di¬ 
rections (not top view) the planet exhibits phases, creating wax¬ 
ing and waning sickle shapes similar to the Moon’s phases. A 
sickle shape gives off less light than the fully illuminated disk. 
Hence the overall (disk-integrated) intensity shows variations 
even though these shapes cannot be resolved ([Green et al.|2003[ 


Dyudina et al.[2005 Snellen et al.[2009 1 . (iii) Eor an inhomoge- 
neous planet surface, the net reflection also changes with time 
because different regions are illuminated at different times as 
a consequence of rotation or orbital motion ( Eord et ar]|2001| l. 
This phase-dependent light curve can, in fact, reveal part of the 
surface map ( [Cowan et al. [2009] [201 3[ l. The study of a planet’s 
surface requires capturing reflected photons, but not spatially re¬ 
solving the planet. 

We advocate making an exceedingly long continuous obser¬ 
vation of one exoplanetary system, as suggested by [Cowan &| 


Agol (20081, for the purpose of obtaining the Eourier transform 


of the net intensity. Ideally, the observation should last several 
orbits and cover a large optical range, without a thermal compo¬ 
nent. It would be possible to use data from the Kepler satellite 
( Borucki|2009| l, or a promising extrasolar system may be found 
by the future PLATO mission ( [Catala|2009| l. The Eourier anal¬ 
ysis of this long-duration campaign has clear advantages. The 
bright glare from the parent star, (i) and (ii) signal from irregular 
time-varying cloud patterns, and (iii) random noise leave a con¬ 
tinuous spectrum and can be filtered out (which is what Eourier 
transforms are generally used for). On the other hand, station¬ 
ary surface patterns from land and oceans on the planet give 
spectral peaks (see the periodogram in ( Palle et ar||2008 1. The 
diurnal (rotational) and annual (orbital) motion generally result 
in a mixed time signal, which could be messy, but (iv) the spec¬ 
trum has nonoverlapping peaks allowing spin-orbit tomography 
(as proposed by Eujii & Kawahara|2012 1. Similarly, in a multi¬ 
planet system, orbital frequencies are usually incommensurable 


1 









































P.M. Visser and F.J. van de Bull: Fourier spectra from exoplanets with polar caps 


Table 1. Planet geometries and viewing configurations 


system 

water- snowball 
world 

• ( 

slushball caps 

C # 

caps no tilt eyeball rings 

« •' C! 

extension 

obliquity 

inclination 

a = 0 

e<f 

<^=1 

e<| 

a = 1 
/3 = Q 

a > /3 
a +/3 < ^ 

a > /3 
a +/3 < J 

6 = 0 

a < 1 

13 = 0 
a + e<^ 

O' < 1 

13 = 0 I3>0 

6=% 

Section 

3.3 

3.5 

3.6 

4.1 

4.3 

4.4 

5 6 


(not in orbital resonance) such that (v) spectral peaks from dif¬ 
ferent planets can also be separated out. Signals originating from 
different effects can be distinguished more easily in the spec¬ 
trum (vi). For example, diffusive (Lambertian) scattering only 
yields even higher harmonics (see [Cowan et al.| 2013| l, which 
could be used to separate it from other effects such as ocean-glint 


(jWilliams &; Gaidos|2008 

Oakley & Cash 2009 Robinson et al. 

2010| Zugger et al. [2010 

'. The signal-to-noise ratio in spectral 


peaks of phase curves is enhanced, (vii) and may become com¬ 
parable to a transit signal (after multiyear integration). 


In this paper, we calculate the Fourier spectra of the light 
intensity from planets with particularly simple, but interesting 
geometries that may arise when the planet is covered only in 
water and ice. The geometry for a spinning planet can be de¬ 
scribed with two circular symmetric polar ice caps, or, when 
a planet is tidally locked, with a circular ocean facing the star. 
These cases are interesting because these planets fall right in 
the habitable zone, where water and ice coexist. The phase light 
curve for these geometries is strongest: Fourier peaks are most 
pronounced for these symmetric cases, and visibility should be 
optimal. We also exploit the property that liquid water predom¬ 
inantly reflects specularly and is otherwise a strong absorber, 
whereas snow is nature’s best diffusive reflector. At the same 
time, a simple symmetric configuration leads to mathematical 
simplicity, with the possibility of deriving closed analytical ex¬ 
pressions for key planetary parameters, such as polar cap exten¬ 
sion a and obliquity fi. 


The structure of this paper is as follows. In Sect.we intro¬ 
duce the general framework of our model; the notation and our 
model assumptions. The influence of fluctuations from the stel¬ 
lar source and of surface-obscuring planetary clouds on the spec¬ 
trum is briefly considered. In Sect.j^we show that the scattering- 
and glint contributions /j and /l to the light curve are generally 
quasiperiodic, and then we demonstrate how the coefficients of 
the periodic spectrum of a banded planet are obtained. We derive 
exact formulas for these coefhcients for the Lambertian phase 
curve and for the standard Fresnel light curve, corresponding to a 
snowball planet and an ocean-covered planet. In Sect, ^we con¬ 
sider planets with polar caps, and derive and study the spectra 
for interesting special cases. The last two sections deal with the 
signals from a tidally-locked eyeball planet (Sect.|^ and planets 
with Saturn-like rings (Sect.|^. These geometries are illustrated 
in Table In the appendices, we discuss how a transit affects 
the spectrum, address the problem of observation gaps, and show 
how the main results of the paper can be extended for noncircular 
planet orbits, by adding corrections due to the (small) ellipticity 
of a Kepler orbit. 


2. Fourier series for iight curves 

The central star is placed at the origin of our coordinate system 
and the planet, with position r{t), moves in a circular orbit with 
time t as 

r(f) = r(f)r = ix{t) + jy(t) - (i cos tot + j sin tot)r. (1) 

We denote with to the (angular) orbital frequency: the orbital pe¬ 
riod is 2.71/w, one planet year. Influences from companion moons 
and other planets on r{t) are also neglected (orbit eccentricity 
e 0 is considered in Appendix [B]). The polar axis h of the 
planet lies in the xz plane and has obliquity (3. The observer at 
Earth is located (at a distance L) in the direction d: 

h - ism/3 + kcosf3, (2) 

d - {i cos cf) + j sin (/>) sin 6 H- k cos 0. (3) 

It is natural to choose the obliquity (tilt angle) (3 and the inclina¬ 
tion (polar observation angle) 6 in the interval [0, |]. 


2.1. Configurations and their symmetry 


We assume that the rotation frequency Q. is sufficiently slow 
that the planet is well approximated with a sphere (for planet 
oblateness, see Dyudina et al.| 2005 1 and neglect precession on 
the timescale of observation. Ice sheets and land masses rotate 
about the axis h. In particular, polar caps do not grow and shrink 
significantly with the seasons. There are no effects from oceans 
waves or foam ( Vokroughlicky & Farinella] I995| l. The atmo¬ 
sphere is at least partially transparent, such that seasonal clouds 
patterns contribute little, although irregular cloud patterns can 
lead to interference with the Fourier peaks. We will describe the 
surface with an (ice-) albedo map and an (ocean-) reflectivity 
map (Ml and Ml, respectively). Because the maps rotate about 
h, they have a periodic time dependence with period 27r/Q. The 
annual and diurnal motion may give rise to periodicity in the re¬ 
flected light signal. Of course, the diurnal period on the planet is 
2n/(Q. - to). 

We neglected non-Lambertian scattering and polarization ef¬ 
fects ( Stam et al.|2006[ Zugger et al.|2010 l. Thermal light from 
the planet is assumed to be absent or filtered out (for the ther¬ 
mal light-signal, see [Cowan et al.|2012| l. The star is also not de¬ 
formed by the planet (for ellipsoidal variations, see Welsh et al. 
20101. We further make the standard approximations that the 
planet is much smaller than the central star s ^ S (s and S are 
the respective radii of planet and star), the star is much smaller 
than orbital radius 5 r, and relative motion of Earth (par¬ 
allax) is irrelevant r L. Since is the fraction of star 

light received by the planet, we normalize the light curve /(f) 
conveniently with F(t) = f(t)Is^ /4-r^, with /, F denoting the 
respective observed stellar and planet luminosities. 
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The reflected light curve for one planet is now in general a 
quasiperiodic function of time (derived in Sect.|^, consisting of 
discrete peaks at multiples of the base frequencies u) and Q; 


/(f) = 2 

«=0 w =0 


(4) 


Since / is a linear function of the maps Mi and Ml, and because 
the Fourier transform is linear, the Fourier coefficients are 
also linear functions of the planet surface maps. ICowan et al. 
( 2013| l calculated the contributions from the spherical harmonics 
T”’(t?, (f) to /(f) (with respect to the rotation axis). Hence, there 
is also a linear transformation between their coefficients Cj" of 
the surface maps and the Fourier peaks /J” defined here. The 
spherical wave expansion would be useful in our context; they 
may almost diagonalize the mapping. 

For zero eccentricity, zero inclination, zero obliquity (e = 
0 = / = 0), a banded planet has no time-varying signal. A non- 
banded planet for the same case, however, can have a diurnal 
signal; then m-peaks appear only for n - 0. If the vectors orbit- 
normal k, the periastron i (in case of an ellipse), observation di¬ 
rection d, and polar axis h (in case of a inhomogeneous surface) 
lie in one plane, then the light signal of a banded planet will be 
an even, periodic function with respect to the inferior conjunc¬ 
tion, when the planet is in-between the star and Earth. There are 
n-peaks occurring only for m — 0 and the /„° are real numbers. 




cos ntot. 


(5) 


n=0 


Otherwise, if these vectors break the orbital symmetry so that 
they are not in one plane, the light signal of a banded planet is 
not an even function, and contains cosine as well as sine compo¬ 
nents. 


2.2. Effect of stellar variability and cloud cover 

Because the reflected light signal involves the product of the 
source function /(f) with /(f), fluctuations in stellar intensity 
on orbital timescales disturb the signal. Consider, for exam¬ 
ple, an ideal light curve consisting of one harmonic at fre¬ 
quency not. If we let dl{t) represent the stellar noise, with 
variance Var/, then the measured signal in the time domain 
is F(t) = (/q -I- + f„ cos na)t)s^ 14r^. The noise term 

dI{t)fo interferes with the wanted signal Iof„ cos na>t. One re¬ 
quires /ol/nl » (Variyi^fo to resolve the cosine. The noise for 
solar fluctuations at 10 ~^Hz is due to the l lyr solar cycle, and 
has a level of IO^^/q (Aigrain et al. 2004 1 . For a quiet solar- 
type star, the scattered light curve c an be observed, in principle. 


Lean 


(2000 1 , for a reconstruction of the 


for l/J » 10-^0. See _ 

past 400yr solar irradiance. Of course, the photometry must then 
be sensitive enough to detect the fraction /„//o of the average 
planet luminosity. It is important to overcome the challenge of 
separating the direct starlight from the planet’s reflection - in the 
telescope, with a coronagraph or an occulting star shade - espe¬ 
cially if one wants to remove large stellar background noise too. 
If the starlight cannot be shielded, we require s^\f„\ » lO^^r^. 
This is possible for Jupiter-sized planets near dwarf stars, as was 
done by Snellen et al. ( 2009[ ). 

An astrophysicist, looking for planets around other stars, can 
search for peaks in the long duration intensity Fourier spectrum. 
In practice, the number of observable peaks in Q will be lim¬ 
ited by stellar noise and observation time (and detector sensitiv¬ 
ity). We now estimate the noise-requirement for observing the 


peaks. For a finite observation time T, one may use the standard 
truncated Fourier spectrum Fj{v) - J F(t)e^‘'"dt; the in¬ 
tegration of t is over [0, T]. For our periodic light curve Q this 
gives 

Ft(v) = ^ -noj- mQ) + (CyiAv + noj + mD.)]. 

nm 

Hence, every peak becomes convoluted with the (truncated) 
spectrum /^(v) of the source /(f). Noise at the frequency v 
from the stellar source is quantified by the power spectrum 
S t{v) - |/r(v)p. For T — > oo, this function has a limit for 
all V 0, with the lower frequencies having the slowest con¬ 
vergence. At V = Q, Sj has a peak of width l/T and height 
S t(0) = Iq^T. Thus, the central peak in It{v) has strength 
To observe the peak /„ at v = nai in F(v), it is therefore re¬ 
quired that /o7’'^^|/„| » SAncdy^^fa (peak maximum above 
noise) and that a>T » 1 (peak sufficiently narrow). We can now 
estimate the noise level at v = ncj, using the property that the 
integrated power spectrum equals the variance Var I. For a fre¬ 
quency interval of size Av > w around the peak at nto, one ob¬ 
tains Warllcj > Var//Av > Srina)). Hence, if a>T > 1, the con¬ 
dition for the spectral noise may be more easily fulfilled than the 
corresponding variance requirement for a time series. We need 
roughly (coTy^^\f„\ » lO^Vo- Hence, nine planet revolutions 
would enhance the signal with a factor 3 over one round trip. 

A similar analysis can be used to estimate the effect of irregu¬ 
lar cloud patterns over a planet surface with, for example, (Titan¬ 
like) glinting lakes. If clouds obscure the surface 90% of the time 
and expose surface lakes 10% of the time, the effect is that the 
underlying periodic glint signal is multiplied with a fluctuating 
function with average /q = .1 and variance Var/ = .1 - .1^ = .09. 
The peaks appear for (a)r)'^^|/„| » (Var/)'^^/o//o = 3/o. This 
requires over 60 periods. Of course, the clouds have their own 
(Lambertian) average contribution to the light curve. 

If sharp distinct peaks are resolved in F(v), the observer may 
obtain a base spectral structure with harmonics /i°, /i®, / 3 ° cor¬ 
responding to the planet’s annual motion. When side structures 
/I'j /a*, etc. of similar pattern can be resolved, this could 
indicate the rotation of a nonbanded planet (the planet has de¬ 
tectable spots). Ideally, the frequencies nat+niQ. have incommen¬ 
surable base frequencies w, Q, because then none of the peaks 
interfere. The Fourier analysis of the reflection signal from exo¬ 
planets will obviously be spoiled if the star variability is strongly 
periodic. Likewise, surface mapping becomes impossible if the 
clouds patterns have diurnal periodicity themselves. 


3. Scattering and reflecting surfaces 

We describe the planet surface with two maps, Mj and Ml, for 
the high-albedo (ice, snow) regions and for oceans (liquid). We 
are interested in high-contrast situations where Mi equals one 
on ice and zero elsewhere (expressions in the next section al¬ 
low for an albedo factor). The ice scattering is assumed to be 
Lambertian, isotropically in all directions away from the sur¬ 
face. Ideally, the liquid map Ml also equals one or zero and in¬ 
dicates whether there is open ocean or not. (Reduced visibility 
from cloud cover could be accounted for with a reduction fac¬ 
tor.) A planet with an ice and liquid surface would give a phase 
light curve that is the sum of the two independent contributions 

/(f)=/l(0+/L(f), 

where I is for ice (methane or water) and L is for liquid (methane 
or water). 
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A light ray from the central star that hits the planet sur¬ 
face has the (unit) direction vector r and can be scattered in 
the direction b of an observer on Earth. With s - ss, we de¬ 
note surface vectors (originating from the planet center). Only 
the surface segment that is the intersection of the illumination 
S' = {s e iSIs • r < 0) (day-side hemisphere) and the visibility 
S" = {s 6 iSIs • o > 0) (visible to the observer) contributes to 
/. The phase light curve, therefore, is a surface integral over the 
spherical lune S' n S" = P. For a diffusive scattering surface, 
the light curve contains two inner products, as in 

The explicit time-dependence in the albedo map Mi{s, t) ac¬ 
counts for the rotation. The intensity is properly normalized: 
using • b) - n, • s) - n, the overall 

scattering for a planet with unity albedo M = 1 is precisely 
d^b flit) - An. The surface S is parametrized with the vector 
function sijj., v) and the surface element in the above integrals 
needs to be calculated with the Jacobian 

X ^|ii//dv. (7) 

We use Miis, t) - Miifi, v, t) for rotating maps (ice, land, ocean, 
or otherwise) in non-rotating spherical coordinates (jU, v); see 
Table for the parameters used. To find the corresponding co¬ 
ordinates {■&, if) on the hxed map Mi(?7, if) in the rotating frame 
with respect to the axis h as function of (ju, v, t), one may use the 
following relation: s • (j x h + ij) = 

^ ^ ^^ p _ ^ sinyS -I- i sin v sin fi. 


3.1. Fourier spectrum from Lambertian surface maps 

We are now ready to derive a general expression relating spectral 
peaks of the signal to a surface albedo map Mj. For a unit 
vector s on the unit sphere in cylindrical coordinates {p, v, z), we 
substitute p - Vl - so that the planet surface is parametrized 
with (v,z), as 


Table 2. Parameters used in modeling 


symbol 

quantity 

t 

time 

2t 

transit duration, or gap duration 

T 

observation duration 

m 

star luminosity, power output 

hiv) 

Fourier transform of star luminosity 

k 

star luminosity average 

Siv) 

power spectrum of /(f) 

L 

distance from Earth 

S 

star radius 

s 

planet radius 

s = ss 

planet surface vector 

r 

planet-star distance 

r = rr 

planet position vector 

a 

orbit semimajor axis 

€ 

orbit eccentricity 

(jJ 

orbital angular frequency, mean motion 

(jjt 

orbital phase angle, mean anomaly 

y 

frequency variable; true anomaly (in App.[B]l 

A 

eccentric anomaly 

n 

rotational/spin angular frequency 

h 

rotation/spin axis 

p 

obliquity, tilt angle of rotation axis 

a 

polar cap extension angle from north pole 

b, c 

inner ring radius, outer ring radius 


ring transmitivity, ring- or ocean reflectivity 

Ml 

ice cap/continent (albedo) map 

Ml 

liquid lake/ocean surface map 

(x,y,z) 

various rectangular coordinates 

(p, V, z) 

cylindrical coordinates in inertial frame i, j, k 

(s,p,v) 

radius, polar and azimuth angle in inertial frame 

(s, 0, ip) 

radius, polar and azimuth angle in rotating frame 

e 

inclination, (polar) observation angle 

0 = 0 

observation is face-on (orbital plane) 

0=f 

observation is edge-on (orbital plane) 

<!> 

azimuthal observation angle 

(jjt = (j) 

phase at inferior conjunction 

Fit) 

planet luminosity 

Fjiy) 

Fourier transform of planet luminosity 

m 

light curve, normalized to 1 for a mirror ball 

fni 

Fourier coefficient for light curve 

g(v, t) 

pattern function, planet map integrated over z 

p/n 

on 

Fourier coefficient for pattern function 


s 


(i cos V + j sin v 


) yjl -z^ + kz. 


It follows that 

-s • r - - Vl - z^ cos(a;f - v), 
s • b - Vl - z^ sin 6 cos(v - f) + z cos 0. 


The circle that borders the visible region is found for s • d = 0, 
hence the limb has the coordinate 


Zlimb(t^) “ 


- cos(v - (f>) 
-y/co?0~+cos^(v^^^ 


( 8 ) 


Using the inner products and the surface element d^s - dzdv, 
the double integral in (|^ is expressed as the following repeated 
integral: 


(jJt- j 

4 r e(v, f) 

“ J dv| cos(wf - 

wt-f- 


where the pattern function g(v, t) contains the z integration over 
the interval between the limb point ([^ and the point with z = 1, 
i.e., 

I 

- J" dz(s • b)^/l - z^Mj(s, t). (9) 

Since the rotation axis n and the observation direction b are con¬ 
stant, g(v, t) characterizes the scattering pattern from the visible 
regions. When v and t are considered independent variables, then 
Miis(v,z), t) is periodic in v (because s is periodic in v) and in 
Qt (modulo 2n) and therefore 

OO OO 

n=0 m=0 


When we make the substitution v - a>t + v 

(i)t— 2 ~ 2 

J" dv I cos(ajt - y)| J" dv | cos v| 


cos nv, 


3n 

2 


4 
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and use this relation in the remaining integral to get 






«=0 m =0 



dv I cos v| cos nv. 


3n 

2 


We hnd that the Fourier coefficients of fi are each related to the 
corresponding Fourier coefficients of g by 


rm _ m 

Jo - ^So > 


rtn _ _ m r, 

- Si’ J2 


-(-1)" 4 


2n 


4n^ - 1 TT 


_ m fin _ Q 

„62n’ Jln+l - 


( 11 ) 

The expressions (fTTli are key results of this paper. The Fourier 
coefficients g™ ofthe pattern function g(v, t) found from the 
surface map Mi(s, t) immediately gives us the spectrum of the 
quasiperiodic light curve of a planet. 

There are several reasons why the linear transformation from 
Mi to /™ is not invertible. Firstly, the region with polar an¬ 
gles 7? > I H- arccos(M • b) below the southern circle (i.e., 
(« • b)^ + {h • < 1) is permanently inaccessible to the ob¬ 

server. Thus any map that is zero outside and arbitrary inside 
this circle gives no signal. (These maps are in the null space of 
the transformation, see Cowan et al.||2013l l Secondly, consider 
the case that Mi is only a function of z, not of v and t. This 
is the case for a banded planet with zero obliquity (see |4.3[ ). If 
a banded map is also zero below the northern circle z = sin 9, 
it can be seen that the integral in expression ([^ for g{v) effec¬ 
tively runs over z-values from sin0 to 1. The dependence on v 
comes from the factor (s • b) only, which yields a single cosine. 
Therefore, the spectrum has just two peaks, with amplitudes /o° 
and /i®, and most of the information from the banded structure 
is lost. Thirdly, if a banded planet with zero obliquity is also ob¬ 
served edge-on, the inclination is 0 = | and zumb = ±1. Then the 
z integral vanishes for -i-l and gives the same overall albedo con¬ 
stant for -1. This implies that the light curve is indistinguishable 
from that of a homogeneous planet. Finally, if the rotation period 
and orbital period are commensurable, which is when Q./a> is a 
rational number, then different combinations of the coefficients 
have the same harmonic frequency and thus interfere. This is 
the case in particular if the ratio is not very large (or very small). 
Different maps can give the same signal, and detailed structure 
is lost. 


3.2. Fourier coefficients for a banded pianet 


The surface pattern of a rotating planet is stationary only if the 
pattern is banded; that is, if the rotation axis is also an axis of 
symmetry for the pattern. In other words, the planet has no spots 
(like Jupiter or Neptune). The function Mi{s, t) is invariant under 
rotation about h and constant in time; then the planet’s rotation 
is unobservable. It follows from (9p0i that g’” - 0 for all m 
except m = 0. In what follows, we only consider banded planets 
(and tidally locked planets, with Q = to), and, therefore, drop the 
index m from all following expressions. 

If we assume the further (reflection) symmetry Mi(—s) - 
Miis), we can simplify the expression for the coefficients g„. Let 
us consider first gi. We make use of the fact that eis oppo¬ 
site for values of v that are ji apart. Points with z > 0, which are 
invisible, are opposite to points with z < 0, which are visible. 
This allows the z integral (j^ for g(v) to be split into a sum of 
integrals for two hemispheres, and to condense into one p inte¬ 
gral over the north hemisphere. Using pdz = -{p^lz)dp, we may 


write 

= ^ fdve-'' fdp<i^«,(S). (12) 

"i l! VTV 

For even coefficients, the phase-factor e^^'"'' for opposite points 
is the same, hence 




dp 


b\p^ 


VT 


Mi(s), n>l. (13) 


The dependence on limb position in (|^ has been removed in 
P2p3|l. Note the distinct absolute-value of the inner product in 


(13 I is absent in ( [T^ . The assumption of reflection symmetry 
must be independently validated before these expressions may 
be used. 


3.3. Fourier coefficients fora homogenous pianet 

Consider the case of a planet with a homogeneous scattering sur¬ 
face. The orientation of the rotation axis is irrelevant and the 
symmetry about the z axis is only broken by the fact that the 
observer is viewing from an inclined direction, and we may set 
(^ = 0 in 0- We calculate the light curve /(f) and its Fourier 
coefficients /„. It turns out that, on the one hand, the time-signal 
/(f) is most conveniently evaluated using a new coordinate sys¬ 
tem that exploits the symmetry about b, but that the peak inten¬ 
sities /„, on the other hand, are most directly derived from the 
function g(v). We thus write the surface vector s in spherical co¬ 
ordinates {s,p, v) with respect to an orthogonal basis hxed by o 
and r as 

(r X d) X d rxb 

s - -cos V sin u -H o sin V sin u H-cos u. 

\r X o| |r X o\ 

We require for (|^ the inner products 

-s • r - (^ Vl - (r • b)^ cos v - r • d sin v) sinju, 
s • b = sin vsinp, 

and the surface element (j^ is d^S = sin p dp dv. Let if be the 
scattering angle, between r and o, so that 

cos ^ = r • d = sin 0 cos cot. 


The condition S • d > 0 implies that p lies in the interval [0, n] 
and f • s <0 implies that v lies in the interval [0, ip\. With these 
definitions, and by setting Mi - I, formula (|^ reduces to a well- 
known form ( |van Hulst|1980| l 

TT tp 

f 4 r r 8sin(,o-8(,ocos(,o 

/i = - I dyu I dv sinyu(-r»s')(s'»o) = - — -. 

0 0 

Completely expanded, this yields 



sin^ 9 cos^ wf - sin 0 cos tot cos ' (sin 0 cos ut) . 


Note that, for top view 0 = 0, the signal is constant. As pre¬ 
viously indicated, for the Fourier components, we substitute 
Ml = 1, / = 0 and (j^ in the pattern function (j^ and obtain 


g(v) - -i Vl - sin^ 0sin^ v - 


cos^ 0 


2 Vl - sin^ 0 sin^ 


-I- sin 0 cos V . 


5 



























P.M. Visser and F.J. van de Bull: Fourier spectra from exoplanets with polar caps 


There are two steps to obtain the Fourier coefficient from this 
function. First, calculate the Taylor series in the powers = 
sin^^ 6, and, then, expand 


sin^^ V : 




2k 

^2FT 


cos 2nv. 


Collecting the coefficients in front of the cosines leads to the 
result 


^0 


SE-4Kcos^e 

37T 


81 


4 sin 

— z —, 82 n+l - 0. 


The functions K — K(sm0) and E - ^(sinfl) are the complete 
elliptic integrals of the first and second kind; the inclination 6 
has the role of the so-called modular angle. The higher-order 
coefficients lead to the hypergeometric function 2 F\{a, Zt; c, x): 


g2n - 


8sin^"0/i\ 
3 22« \n] 


2 F\{n+ 2 ,n - j',2n+ 1; sin^ 0) 


4sin^"0cos2fl 



2 F 1 {n+ \,n+ 2 ; 2 nH- 1; sin^ 6). 


The hypergeometric functions in g 2 n can be expanded into K and 
E . With ( [TT] ), the peak values /„ follow. We obtain for the first 
four coefficients 


/o 

/2 


32E-l6Kco^^6 


/i -- 


4sin0 


(14) 


32K{2 - 3 sin^ 0) cos^ 6 - 64E cos 20 
21 sin^ 0 


/3 = 0. 


If one observes the orbital plane edge on, then 0 - \ and 
g>(t) - \<jjt\ for -n < cot < n. The standard Lambertian phase 
curve for in-plane observation of a planet or moon is 


Mt) = 


8 sin Ifatfl - 8|wf| cos cot 
3n 


—n < cot < n. 


(15) 


with Fourier coefficients. 


32 -4 64 1 

37r2’ 3’ 37r2 (4n2-l)2' 


(16) 


The graph is plotted in Fig. The coefficients ( [T4| ), ( [Thl l drop 
to zero with the fourth power in n, such that /4 is already much 
smaller than |/i|. Fig. shows the ratios of coefficients g„ as 
function of observation angle 0. If the planet radius s and albedo 
were determined with other methods, then the first coefficient fi 
could be used to estimate or verify obliquity. Otherwise, the rel¬ 
ative values of the peaks in the spectrum must be used to obtain 
orbital parameters. For example, the ratio of the first two peaks 
/ 2//1 = -4g2/3ngi could be used to find 0. For this purpose, in 
Fig. I^we plotted the ratios as a function of 0. The overall com¬ 
plex phase <p of the peaks, which was set to zero for the sake of 
convenience in our calculations, is just the orbital phase with re¬ 
spect to the inferior conjunction (in-between Earth and the star). 
The absence of /s, /s would indicate that the planet is a diffusive 
scatterer without specular glint. The condition \ f„\ » 10^^ fo, de¬ 
rived in Sec. |2.2| valid for a sun-type star observed for one year, 
implies that /g is swamped in the stellar noise. 



Fig. 1. Standard light curve (|15| for a homogeneous 


(Lambertian) planet; normalized intensity / versus orbital phase 
cot, for edge-on observation. The inset shows, to scale, the coef¬ 
ficients g„ which are proportional to /„ according to ([TT]i. 


g2n 



Fig. 2. Harmonic ratios g 2 /gi (blue), g 4 /gi (purple) and g(,/gi 
(orange) for a snowball planet, versus inclination angle 0. These 
curves can be used to determine the inclination 0 or to conclude 
that the planet is, instead, inhomogeneous. 


3.4. Glint spot on a specular reflecting surface 

Consider a reflecting sphere with radius s, placed at the origin, 
which is illuminated by a point source on the z-axis at z = -too. 
Imagine a long cylinder of parallel light rays around the z axis 
with radius p < s and cross-section A = np^. The rays inside this 
cylinder hit the north pole of our sphere from angle zero up to an¬ 
gle 0. Since s is the radius of the sphere, this implies p = s sin 0. 
At the border, the angle of incidence is also 0, so the cylindrical 
beam is reflected into a segment of the hemisphere with limit¬ 
ing angle 20. The total solid angle of the reflected light is there¬ 
fore given by the integral fg dt^sint? = 4AIs^. Because 

this solid angle is proportional to A, a perfectly reflecting sphere 
scatters light isotropically in all directions. Since a finite source 
(at large distance) may be composed of many point sources, a 
reflecting sphere also scatters the light of a finite source isotrop- 
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ically. The reflection of a metallic planet would, like a garden 
ball, be equally bright from all directions, and there would be no 
variation in intensity during its yearly orbit. 

If the light source is a star with radius S and our sphere is a 
planet with radius s, the light rays that are incident on the smaller 
planet are conhned to the cone rr + rdr. Here dr is orthogonal to 
r and its length is |dr| - S fr. Since the distance to the observer 
is very large, the reflected cone is practically collapsed into a 
line with direction o - Ld. The points where the light rays are 
reflected on the surface, form a small ellipse around the central 
point of reflection 


, o-r o-r 

s — -=-. 

|d-r| 2 COS)? 

Here )? is the angle of incidence with respect to the surface nor¬ 
mal, with cos § - -s • r > 0. This ellipse is now given by 

d-r-dr b- r (s • dr)s - dr 

s H- ds = -=-1-. 

|d - r - dr| |d - r| 2(-s • r) 


The vector ds points from the center to the edge of the ellipse 
and is orthogonal to s. The ellipse has semiaxes of S/2r and 
S/2rcos§ and an area of d^s = 7r5^/4r^ cos)?. Ocean waves, 
however, increase the size of this glint spot. As discussed by 
Williams & Gaidos (2008|l, the reflection signal does not signif¬ 


icantly change as long as the ocean is wider than the spot. In the 
limit S/r —> 0, the image of the star becomes a point (except 
for s • f - 0, which is the case of grazing reflection). The con¬ 
clusion is that, in the regime S r ^ L, the reflection signal 
equals the surface reflectivity at the point of reflection and thus 
probes the surface along a curve in the band with polar values 
^ - jO < /j. < f + 50 . Close to edge-on observation 0 < |, the 
entire northern hemisphere is probed because, near inferior con¬ 
junction, the glint spot quickly moves from the yu = | to = 0. 
If we call ??()?) the ocean reflectivity for incidence angle ■&, the 
signal is 


Here and Ml(J, t) is now the (rotating) ocean map. When taking 
a fixed vector s on the surface, this map has the Fourier series 
Ml(s, t) - Xm=o Re(e''”^'M'"(S)). The glint point (in the nonro¬ 
tating system) varies periodically over one year as s(t + 2nlco) - 
s{t). When we take this periodicity into account, the coefficients 
for the rotation map are so that the net 

function f) is a quasiperiodic function of t. The reflectiv¬ 

ity ?? is a function of time via §{t), and also varies also annually. 
Therefore, the light curve for ocean glint, is a product of a pe¬ 
riodic with a quasiperiodic function, which is quasiperiodic and 
has the form Q, just like the light curve for a Lambertian map. 


3.5. Fourier coefficients for ocean giint 



Fig. 3. Standard Fresnel light curve ( |18| from ocean glint for 
edge-on view (in blue; the Lambertian scattering signal is shown 
for comparison in gray.) The inset shows the spectrum. In the 
visuals, the glint is yellow. 


Taking n = 4/3 for the index of reflection for water (or liquid 
methane), and assuming that the stellar light is unpolarized (the 
relative s- and p-intensities are ^ and j) the net signal from a 
smooth water surface is 

^_??s+‘Rp_ 49(9-llcos2)?-t 18 cos4)?) 

^ (21 -H 75 cos2)? H- 25 cos)? V7 -t 9 cos^ 

(17) 

For a planet in orbit, with position vector Q, observed from 
direction ([^, the cosine of the incidence angle varies with time 
as 


cos )?(f) = -r»s-s»d = 


^ ■ 


This must be substituted in (17 1 to get the light curve for a water 
world (with ocean map Ml(S, t) = 1). The signal is significant if 
the observer is in the plane of the orbit, for {6,(p) = (f,0), and 

cos )?(f) = Vl - cos lotj V2 = I sin ^cot |. For this case, we obtain 
the standard Fresnel light curve, which is 


98(16 - 7coscuf H-9cos^ wf) 

/l(0 = - - - 2 - 

(l 17 - 75 cos ut + 25 V46 - 18 cos (jjt\ sin 

(18) 

The glint spikes at f = 0, for grazing light rays, when r and o are 
in the same direction; see Fig. The width extends to about a 
quarter of a year. The numerical values of the Fourier peaks are 


In this section, we derive a compact expression for the ocean 
glint, the ideal Fresnel curve, and calculate its Fourier peaks. 
The reflection by a surface of water depends on the angle of 
incidence )?. The reflection coefficients for s-polarized light and 
p-polarized light are given by the Fresnel equations ( Hecht|2002|l 


/o = 0.12032, (19) 

/i= 0.181, /2= 0.141, /3= 0.102, 

/4 = 0.075, /5 = 0.055, fe = 0.042, and 

f„ =» 50^3 if n —» oo. 


W) = ( 

!^p(t?) = ( 


Vn^ - sin^)? - cos)? 


yjrfi - sin^)? + cos t? 
n^cosi?- 

rfi cos § + ^JrP- - sin^ ?? 


)■ 


Figure 4 shows the ratios of peaks that do not occur for a 
Lambertian planet. Observation of these values would strongly 
indicate that the planet is covered with oceans, composed of wa¬ 
ter or methane. The coefficients are positive and decay with an 
inverse square as /„ ~ which is characteristic for func¬ 
tions that are differentiable except at one point (in this case at 


7 
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t - 0). Many peaks can come out of the noise, in theory, for 
\f„\ » 10“^/o. However, the actual sharpness of the f = 0 peak, 
will limit the validity of our calculated numerical values of /„. 


3.6. Slushball Earth 


An interesting situation occurs for a planet that has an ocean cir¬ 
cumventing the equator, and is covered with ice everywhere else. 
For Earth, this is called slushball Earth ([Hoffman et al.||1998[ 


Fairchild & Kennedy|2007 Macdonald et al.|2010 1 ; it is the hy¬ 
potheses that during the periods of global iceage, almost all of 
the oceans were covered with ice sheets, except for a narrow 
band at the equator. For an edge-on view, the signal has a max¬ 
imal intensity; there is the scattering ice with nearly unit albedo 
and there is a maximal glint effect. The planet can at the same 
time be found in a transiting orbit, i.e. eclipsing the host star. The 
light curve is found with the summation of ( [T5] l, (18 1 , 

Fourier coefficient is the sum of the expressions (|16|l, (|19|l, ( |A.2| i 
that follow. 


4. Polar caps 


fln+l 



Fig. 4. Harmonic ratios fs/f^ (blue), (purple), and /g/Zs 
(orange) for a planet with permanent ocean-glint. The glint spot 
of a water world with polar caps is continuously visible for 9 > 
2a+ 213- f. 


We describe the situation with polar caps as 


Miis) - H{\s • h\- cos a), Mi^(s) = 1 - Mi(s), (20) 


using the Heaviside step function N, assuming that the caps are 
constant in size. Seasonal growth and decay of the caps is as¬ 
sumed to be small, which is reasonable for larger planets or small 
obliquity. Hence, Mi - Q for - cos a < s • h < cos a, and unity 
otherwise. (An albedo factor and a cloud reduction factor could 
be included.) The border of the north polar icecap is found where 


cos a 


- s • h - psm/3 cos v + -^1 -cos/?. (21) 


We assume /? < a, so that point sk lies on the north polar cap and 
we assume a + f3 < f, so that the tropical latitudes are ice free. 
The radial coordinate p for the polar cap edge cap, as a function 
of V, denoted pn, is then the positive solution of the quadratic 
equation (211. One may write this equation and this solution as 



2 cos^jS “ cos^ a -I- 2pcos asinjScos V 
cos^ jS -I- sin^ [3 cos^ v 

cos a sin (3 cos v -i- cos /? -^sin^ a - sin^yS sin^ v 


( 22 ) 


Pn(v) 


cos^ P + sin^ p cos^ v 


Fig. 5. Light curve of planet with polar caps for edge-on obser¬ 
vation, (0,0) = (f>0)> with and without glint (blue and black 
curve). The caps extend to a = | and the axial tilt is toward 
the observer with P - Inset shows the spectra, also with and 
without glint (blue and black). 


The border of the south-polar cap is found at Ps(v) = pn(v -i- tt). 


4.1. Case I: edge-on view 

Consider the case that the orbital plane is observed edge on, 
when 0 = I. This important situation arises for extrasolar planets 
that were discovered with the transit method. During an eclipse, 
the planet passes between its parent star and Earth, so we ob¬ 
serve the orbital plane edge on. The limb is no longer given 
by (IHl), but we hnd that v = 0 ± f • Between these values, for 
(f) -j < V < (j) + j, the pattern function is 


g(y) _ 1 

(4 2+p2 

.2 2 +A / 

2 ' 

is 3 > 

1 Pn 3 V 


yjl -p2jcos(v-0). (23) 


Substitution of (20 1 in (T2p3]l yields the Eourier coefficients 


TT 

gi^^J dv ^ -y/l - p2 j cos(v - 0), 

-7T 

n 

gin ^1 -p^jl cos(v - 0 
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Fig. 6. Light curve of planet with polar caps for edge-on obser¬ 
vation, {0,(p) - (j, ^), as in Fig.|^ except with the axis orthogo¬ 
nal to the observation direction. 


Here we need to substitute the cap-edge function (22 1 . For small 
y8, the coefficients are 


Fig. 7. Light curve for planet with polar caps for inclined obser¬ 
vation 0 — The caps extend to a = and without axial tilt 
(3-0. The north polar cap is completely in view, leading to a 
pure cosine. The glint is interrupted by the cap, but the effect is 
small. The gray curve is the Lambertian signal for this configu¬ 
ration. 


o -1^2 -H cos a . 4 a 

£i = 8 e ^ -sin — 

SI 3 2 

7,74 2 -H COS a . 4 

g2 = 32e-2'(»—^-sin^ 

97T 

15 -I- 5 e- 2 i'A _ 


+ 20^ 


IStt 


-I- 0 sin a sin 2a cos 0, 

a 

2 

sin a sin 2a. 


If transits occur, the large signal ( |A.2[ ) from the stellar light dip 
has to be carefully subtracted first. Figures and show light 
curves for an edge-on view in the respective cases of 0 = 0 and 


4.2. Case II: one polar cap in full view 

When the limb does not cross the icecaps, one cap is always 
completely visible and the other cap is entirely obscured. This 
allows us to go from z to p in equation (|^ and obtain, with sub¬ 
stitution, ( |20l i that 


g(y) 

2 


Pn 

/dp( 


VTV ' 


(24) 


(2 2+pI 


pI 


■\Jl -pnj cos(v - (/i) sin 6 -I- cos 0. 


For small fi the first coefficients become 

16-1-8 cos O' . 4 a . . . 

= e ^---sin — sin 0 -I- p sin a sin 2a cos <p cos 0, 


g 2 — ^ sin^ a sin 0. 


that has full view onto the polar icecap. Then - p 
and ( |9|24| i, with (p -0, become 

£(v) (2 2 - 1 - sin^ a \ sin^ a 

—— = -T-cos a cos V sin 0 H — 


s = sin a. 


cos 0. 


There are only two nonzero Fourier coefficients: gQ and ^i. The 
relevant nonzero Fourier coefficient in the signal, /i = is 
reduced from -4/3 by the factor 


a 


(4-1-2 cos a) sin — sin 0. 


Figure shows an example of this case. A short calculation 
shows that the glint spot is on an open ocean for 


cos ut < 


^(1 

sin0'' 


cos^ 0 
2cos2 a 


The case of zero obliquity with the observer inside the orbital 
plane gives the standard Lambertian curve (151 with harmonics 
( [T6| |, as in Fig. except that Eqs. ( [T5p6 1 must be multiplied 
with the cap size dependent reduction factor 


(4-1-2 cos a) sin"^ 


a 

2 ' 


Of course, for an edge-on observation of a water world with a 
polar cap, the ocean glint is maximal. This effect turns up in 
the odd harmonic amplitudes, f 2 n+\, of the spectrum. When the 
reflection signal is subtracted from the spectrum, the Lambertian 
contribution remains, and the extension angle a of the polar caps 
can be determined from that contribution. 


4.4. Case IV: face-on view 


4.3. Case III: zero obliquity 

If the rotation axis is vertical, the case p - 0, the symmetry 
about the z axis is only broken in the observation direction ([^. 
Consider the case that the planet is observed from a direction 


We consider, as a last concrete example, the situation that the 
observer has a top view onto the orbital plane, with the entire 
northern polar cap also in view. For this case, 0 = 0, and we 
require that a-\-p < |. The pattern function (24 1 for this situation 
in terms of the cap edge (22 1 is g(v) - 2pJJ/3. To find the light 
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curve itself, the Cartesian coordinates (x, y) are most convenient 
as parameters for the northern hemisphere, where 


s - -rx + (r X k)y + kyjl - x^ —y'^. 

The X and y coordinates are now the corotating coordinates, so 
that the day-night terminator circle found for s •r = 0 is at x = 0. 
The surface element 0 is 

, 1 

d S = — — dxdv. 

Vl - -y2 

and the inner products are -s»r-x, s»d — yjl - x^ - y^, 
which implies that (-s • r){s • d)d^s - xdxdy. The border of the 
polar cap is located at s • n = 

cos a - -^1 - x^ -y^ cos + {y^inojt - x cos ojt) sin (3. 

The quadratic form for the projection of the ellipse onto the xy 
plane is 

(y^ H- ;ii:^ - 1) cos^/1 + (cos a - y sin/1 sin cut H- xsin/3cos = 0. 


If we solve y as a function of x, we find solutions y+(x) and y_(x). 
The extreme values for x, where these two solutions connect, are 


x± - ± sin a -^1 - sin^jScos^ cut - cos a sin /I cos wf. 
The surface integral in (j^ can now be calculated using 


I 


dx (y+ - y_)x : 


4 cos p sim a 
(x+ - x_)2 


J" dx V-^+ - xyfx 


The light curve is therefore given by 

12x^ — 8x+x_ + 12x^ /- , 

flit) = - 3 ^(x ^2 - cos/3 sin^ a (25) 

/ 2 X+ H- X_ , \ ^ . T 

-H - arctan —+ 1 (x+ + x_) cos p sin a. 

2 yl\xfxl\ ' 

The graph is plotted in Fig. for polar caps that extend 45 de¬ 
grees and for obliquity of 30 degrees . 


5. Tidally locked planet 

When a planet is tidally locked, the same side faces the host star 
at all times. The substellar point s = -fs is constantly directed 
toward the star and the point s - rs on the other side is always 
in the dark. Apart from moving cloud contributions, the surface 
pattern will be stationary in the corotating frame. A reflection 
signal from the planet is time dependent because an observer 
sees a rotating surface pattern with a constant illumination, in¬ 
stead of a time-varying illumination, as before. Since this is the 
special case Q. - o), the general equations of Sec.j^are valid, but 
now all coefficients contribute to the same noi-peak in Q, 
for all m. We proceed, however, by defining a new pattern func¬ 
tion. For this situation, we introduce rotating coordinates (v, z) 
defined with 

s - i-r cos V - X r sin v) + kz. 

In this coordinate frame, the surface maps Mj and Ml are only 
functions of (v, z), and they are not functions of time. The planet 
will now have a rotationally symmetric map with respect to the 



Fig. 8. Polar-cap light curve ( |25| for top view; inclination 0 = 0. 
Cap extension is a = |, obliquity is /3 = |. Graph and spec¬ 
trum are clearly different from the situation in Fig. The glint 
is now constant and can be only be inferred from the caps inter¬ 
ruptions. The gray line is the constant intensity from a homoge¬ 
neous planet for top view. 


time-varying f axis, but not with respect to any fixed axis. One 
expects that for 0 4^ 0 and 0 4^ f, higher-order odd harmonics 
from both Mj and Ml are present in the spectrum. For edge-on 
observation b - i, however, a similar symmetry as for a banded 
planet with a fixed rotation axis arises. We now have 

-S • r - Vl - cos V, s • b - - ^jl - cos(wf + v). 

Equation (j^ may therefore be expressed in terms of a pattern 
function g for a locked planet, according to 


4 r g(y) 

/(f)=-J dy|cos(cuf + v)|^. 


g(y) 

2 


1 




cosv I dz(l 
-1 


■ z^)Mi(s), for 


TT 

- < v< 


TT 

2 ’ 


and g(y) is zero for the dark side of the planet, |v| > |. If we 
take the same definition ( [TOl l, then the relation ( [TT| l for m - 0, 
between the harmonics in the signal function / and the pattern 
function g is valid. 


5. 1. Eyeball worlds 


Consider now a tidally locked planet with open water. At the 
star-facing side, a circular ocean exists with an extension an¬ 
gle a ( PieiTehumbertj 20TT] l. The rest of the planet is covered 
in ice, as depicted in the vi suals of Fig. [9] Be cause of oceanic 
heat transport, according to ( Hu & Yang|2014| l, the ocean shape 
may deviate from a circle for the relevant case of a locked planet 
near an M-type star. In our model, the edge is now located for 
-s • r = cos a. We therefore assume the following form; 


Mi(S) = H(-s •r- cos a), Ml(S) = 1 - Mi(S). 
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Fig. 9. Light curve ( 26||27 1 for an eyeball planet with circular 
ocean extending over a = njA, with and without glint (blue and 
black), as given in Eq. ( [T8] l. The view is edge-on. The snowball 
light curve (gray) is for comparison. 


fn 



Fig. 10. Peak ratios fjfo (black), / 2//0 (blue), / 3//0 (purple), 
and / 4//0 (orange) for an eyeball planet viewed edge-on plot¬ 
ted as a function of the ocean extension angle a. Ocean glint is 
neglected. 


Fig. 11. Light curves (28i, (30i from Saturn-like rings with 
obliquity /3 - ^ and H - T - \ (so that rings can be ob¬ 
served lit-from-behind) and zero inclination; with and without 
planet-casting shadows (black and gray). Then, shown in the in¬ 
set /i = 0. For orbital phases u)t - ±|, the rings are edge-on 
illuminated and the signal is zero. 


visible, and the phase curve is 

8 sin a)t - Scat cos cat 
flit) =- 


371 


37r 


/ / cos^otn r 

dv| cos(wf -H v)|(2 + -"Vi 

^ COS^V'^ 


2 -cot 

8 sin cat 


37T 
8 t 

“sL'" 
8 


COS ut 3 cot cat \ 

-h cos a arccos- 

sin a tan a ' 


(26) 


cos cat 


® 2 9 2 ^ ■ 

- — V sm wf - cos^ a cos wf - — st 
37r 


37r' 


sin^ cat - cos^ a 


\ 3/2 

) ■ 


For cat > f + a, the entire ocean is visible, and hence we need 
to subtract the above integral over the entire interval [-a, a]. We 
find 


flit) = 


8 sin cat - Scat cos cat 


3n 


• (1 - cos^ a)- 


8 cos cat 


(27) 


We equate g and /. We have for a < |v| < | that g{v) — 
(8/3) cos V and for |v| < a that g(v) - 


For cat > n - 2a, the glint spot is visible and a contribution 
/l, given by ( [T8]l, must be added to the signal /. This time sig¬ 
nal, shown Fig. 1^ has a peculiar wiggle in the interval when the 
ocean is in full view. Clearly, the pattern function g is most use¬ 
ful for evaluation of the spectrum. The leading term in the large 
n expansion is 


4 cos V 



8 4/ cos^crx , -r--— 

- cos V-(2 -I- 3 — 1 Vcos^ V - cos-^ a. 

3 3 ' cos^v'^ 


/„ ai V2sin2a(-1)" cos(2na - if n —» 00 . 

If one wishes to determine a from the Fourier peaks, and verify 
that the signal is from an eyeball planet, the curves in Fig.[T0|can 
be used. 


At V = ±a the pattern function has an infinite respective left- or 
right-derivative, leading to ~ 

For 0 < wf < f - a the ocean is not visible, and / equals the 
snowball signal (flSll. For j — a < cat < j + a the ocean is partly 


6. Saturn-like rings 

Consider a Saturn-like planet, a planet with a large disk-shaped 
ring system. We model the ring system as a homogeneous and 
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the thin disk in the plane with normal vector h as in Q. With 
the inner radius b and the outer radius c, the total surface area is 
nc^ - nb^. We study conhgurations with s <ih» d)b, so that the 
planet does not block a part of the observer’s view on the ring. 
Light scatters off the ring with reflectivity K. For a thin disk, 
light can pass through so that the rings are also visibly back lit 
(see |Porco et al.|2005| l. The light that comes out the other side 
is again scattered diffusive, with transmitivity T. This approach 
requires the disk to be optically thick (for optically thin rings, 
Arnold & SchneidS] 2004|l. We factorize the net intensity 


F{t) for the light scattered off fhe ring according to 


F(f) = 


• nb^ 


Anr^ 


+ T)If(t)(h • d) sin/?, 


with /(f) independent of the ring absorption \-%-T, and also 
with the dependence on observation direction o factored out, as 
well as the sine of the ring tilt angle /?. 

Consider hrst the simplihed case where the planet is much 
smaller than the ring disk, hence s « c, so that the planet’s 
shadow is unimportant. For the interval \u)t\ < | where the illu¬ 
minated face has normal s - -h directed away from the observer 
and only transmitted light is visible with intensity T, whereas 
for the other half orbit, the illuminated face has s - h, with the 
scattered intensity /?. The phase curve for the rings is then 


/i(f) = 


I cos cut I -I- cos ut T 


I cos (jjt I - cos cot 


2 H + T 

T-% 1, 

—--cos (jjt ^— cos a>t\. 

1% + lT 2 


<R 

K + T 
(28) 


The Fourier coefficients are 


1 T-K -(-!)" 2 

/o - /i - 2<R + 2T' “ 4n2- 1 - 0. 

(29) 

The case of opaque rings, for T - 0, actually produces the same 
signal as polar caps in the small-cap size limit a —> 0. When 
0 < a « jS the polar cap is located inside the polar circle of 
latitude. The surface integral (j^ is reduced to a very small circle 
around h so that the factor {s • o) can be considered constant. 
The signal is a cosine function when this cap experiences polar 
day, and the signal is zero during polar night. This situation is 
opposite to the cases from Sect. where we only considered 
a > p. 


6.1. Shadow cylinder 

We discuss, as a final application, the situation where the planet 
also casts a shadow on the rings. In this case the planet itself 
and the shadow cast by the rings onto the planet are neglected, 
compare [Arnold & Schneider| ( |2004| l. This shadow is obviously 
important when the planet is not small compared to the rings. 
Consider the case when the planet is so large that the shadow 
cylinder of the planet always intersects the ring and at any time 
part of the ring is dark. This means that for our calculation the 
obliquity is bounded by 


c sin/? < s < b cos(ff - /). 

We calculate the illuminated effective area using the plane or¬ 
thogonal to r. If we view the area from the position of the cen¬ 
tral star, this area is precisely the region not obscured by the 
planet disk. This geometry appears in Fig. 12 (on the left). We 
hrst imagine that the ring starts at the planet surface, for radii be¬ 
tween s and c. Since the real system has radii between b and c, to 




Fig. 12. Ringed planet, with planet radius s and outer ring radius 
c viewed from the direction -r of the star, dehning / (left side) 
and from the polar north h, dehning the angle (right side). 
The illuminated area consists of A plus C. On the left, A is the 
difference between two half ellipses with areas A + B and Z?; and 
C the projected circle segment C + D minus a circle segment D 
(with respective areas sin/ and (f>s^). 


obtain the required signal we must subtract the same expression 
with c replaced by b from this initial result. Thus consider now 
hrst a ring running from s to c. The cross section of the half of 
this ring that is on the dayside of the planet is (c^ - ^ 2 ) 1 ;^ • h\nl2, 
since the dayside is fully illuminated. The illuminated part of the 
nightside half of the ring consists geometrically of two (identi¬ 
cal) segments of the circle of radius c with angle projected into 
the outer ellipse, minus two circle segments of radius s with an¬ 
gle / (see Fig.[T2|i. Hence, / is the angle between the point where 
the shadow circle of the planet crosses the outer (projected) el¬ 
lipse and the direction normal to h, whereas is the angle be¬ 
tween these points for a face-on view. This yields a cross-section 
c2|r • «|/' - *2^. Substitution of the point on the circle in the 
equation for the outer ellipse gives 

/scos(/i\2 /ssin(/i\2 ^ 

^ c ' ''cr • h' 

It follows that the angle is given by 


(p{c, t) - arctan 


2 2 
C - i 


^2/ C0s2 OJt sin2/ - c2 


The relation between and / is |r • h\ tan <p' - tan <p. Hence, the 
second angle is given by 


4>'{c, f) — arctan 



If we now add the contributions of the two ring halves, we obtain 


. 1 , 1 0 - 01 
flit) = - I cos Wf| -H-r- - - I cos Ut\ 

1 7T 

^2 (pic, t) - 4>{b, t) 


- b^ n sin /? 


where we must remark that these equations are only valid for 
\cot\ > I and T - 0. This is the situation when the rings are 
opaque and only visible half the time. The general case is just a 
linear combination similar to (|28|l. To lowest order in the axial 


12 































P.M. Visser and FJ. van de Bull: Fourier spectra from exoplanets with polar caps 


tilt y6, and for b - s, the result is simply proportional to Equation 
( [28) l, but with an extra geometric factor 

1 s s 

2 n{c^ - s^) c 

When ( p9| ) is multiplied with this same factor, the asymptotic 
behavior of f„ for large n is accurately described for arbitrary 
p. For the case of maximal visibility, and arbitrary /?, which is 
h - b,b - s, and c - si sin^S, the signal from the rings simplifies 

to flit) = 




2 ncos^p 


arctan 


COtyS ' 

I sinotfl, 


I cos wf| 


sin/? 

TTCOS^P 


arctan 


cos/? 

I tan (jjt\' 
(30) 


The graph appears in Fig. 11 The signal is significantly weaker 
because for b - s, the planet shadow is large. 

The spectrum of the signal from a planet with rings is of 
the form Q and is again characterized by the fact that higher 
order odd coefficients are all zero. For T also /] = 0. The 


dominant Fourier coefficient of (30 1 is a relatively short formula. 


h 


1 


blp cos^/? 


^3;7r + TTCOS 2/? ■ 


4 8 log cos /? 

tan p tan^ /? 


•4y6). 


\n+l. 


The Fourier peaks for rings behave asymptotically as (-1)" ‘n 
like (29 1 , whereas for the planets the tail behaves as like 
This behavior allows one to discriminate between the sig¬ 
nal from a thin disk and a sphere. 


7. Conclusions 


When a planetary system is observed continuously over a num¬ 
ber of years, the Fourier transform of the light signal provides 
a natural method of analysis. The immediate advantages are: 
(i) The glare from the host star (in observations for large planets 
near M-dwarfs without an occultor), or some residual glare, can 
be filtered out. Stellar noise at the orbital frequencies swamps 
the Fourier peaks for small (ii) Atmospheric cloud cover 

that partially obscures the surface, is also partially filtered out, 
if clouds mainly contribute an a-periodic component in the sig¬ 
nal. (iii) Well-defined orbital and diurnal frequencies (see Palle 
2008|l can be obtained from the peak-center positions. 


et al. 


(iv) Contributions of different planets can be identified, provided 
they are not in orbital resonance, (v) The asymptotic behavior of 
the peak-amplitudes, phase factors, and even-odd structure al¬ 
lows us to separate out and identify different contributions (re¬ 
flective oceans, land or ice sheets, or rings) in the light signal, 
(vi) For observation time T of N orbits, the signal-to-noise ratio 
in the peaks is enhanced by the root of N, essentially because the 
Fourier transform is an averaging method. 

How can one distinguish between planets with or without 
rings and identify the different configurations we consider here? 
Since none of the specific parameters are known in advance, one 
has to rely on general properties of a measured Fourier spec¬ 
trum. To apply our results, an astronomer should first identify 
the spectral peaks from the intensity spectrum of a star system 
that belong to a single planet. If sidebands to the annual ncj- 
peak pattern (a similar pattern shifted by Q and 2Q, etc.) are not 
found, either the planet is banded and has no longitudinal pat¬ 
tern, or the planet is tidally locked (so that Q. - a>). Then the 
astronomer needs to fix the orbital phase with respect to the in¬ 
ferior conjunction. For a transiting planet, this is at moment of 
transit (peaking). 


A generic Lambertian signal from a 
If we assume a stellar-noise level 


If the coefficients are all positive and decay as /„ ~ this 
may indicate ocean glint, see Eq. ( [T9| . If there are no higher- 
order odd peaks, fin+i = 0, then there is no glint, which implies 
that the point of specular reflection is on land or behind a thick 
atmosphere, see Eq. When the glint spot abruptly moves 
from ocean to a polar cap or back, this generally occurs when 
the glint effect is weak. We therefore expect that the behav¬ 
ior due to such a discontinuous jump will only dominate in the 
far-end regions of the tail that cannot be resolved. If the spectrum 
(after subtracting any observed glint component) has an alternat¬ 
ing even tail behaving as / 2 „ ~ (-l)"^'n“^, the signal is that of 
a scattering disk, either from a flat polar-cap (inside the polar 
circle) or from a ringed system. The value of /i can determine 
if the signal is light scattered by semitransparent rings, see Eq. 
( |28| l. The peaks of an locked planet with a circular ocean decay 
in strength as \f„\ ~ 
planet decays fast as / 2 „ ~ n 
of 10“^, the number of peaks that is expected to be resolved is 
about 100, 30, 10, and 6, for a perfect glint spike, rings, and 
Lambertian scattering ice caps, respectively. If the photometry 
is not good enough to detect the fraction 10“^ of the average 
planet luminosity, of the order i^//4r^, then fewer peaks can be 
found. 

For observation that is not edge on the orbital plane, the spec¬ 
trum obtains an extra factor of sin^" 6. This exponential decay 
will limit the number of peaks that can be resolved, which will 
hamper identification. The effect of orbit eccentricity, see Eq. 
Q is also exponentially weak because /„ ~ (ee/2)". This ef¬ 
fect dominates only for zero inclination, and in case obliquity is 
zero or the planet is homogeneous (or at least the exposed re¬ 
gion is). Any other periodic signal is caused by changing phases 
or rotating patterns on the planet surface. The expressions at the 
end of Sects. 4. IE l^for gi, g 2 are simple enough to see whether 
the four parameters, 6, a,p, for a water world with polar caps 
fit the observation of the first two Fourier peaks. 

If one isolated peak appears in the spectrum, this could indi¬ 
cate a homogeneous but large planet at such a small inclination 
that all other peaks are suppressed. However, if all the reflected 
light comes from a circle around the (north or south) pole that 
happens to be in full view, a pure cosine light curve is precisely 
what will be observed for a planet with zero axial tilt, see Sect. 


4.3 A bright circularly symmetric albedo map around a planet’s 


pole on a further absorbing surface, could indicate an exposed 
circular polar ice cap surrounded by oceans. 

The main disadvantage of the Fourier transform is the need 
for many periods, hence a long continuous observation (in the 
optical range), but we think that it would still be worthwhile to 
single out a promising exoplanetary system as a candidate for 
study of the intensity spectrum. 
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Appendix A: Transits and observation gaps 

When planet transits and reflected light curves are observed in 
one signal, large amplitudes from the periodic transits are added 
on top of the weak reflection peaks. For Jupiters near dwarf stars, 
it may still be possible to separate the two components in the 
signal. 

Observing from direction d, one sees a planet passing in front 
of the star disk for |r x d| < S. The transits are thus found (taking 
ip - 0 for convenience, and r 5, and 6* |) in f-intervals for 

which 

r sin \u)t\ < r sin wt = ''Is^ - cos^ 0. 

Then the measured signal dips by -s^IqIS^ (which is -10“"^/o 
for a twin earth). We are interested in the reflected signal from a 
planet. This is a factor 5 ^/4r^ weaker (10^®) than a transit signal 
and weaker (10“'°) than the star. For a Jupiter close to 

a dwarf star, these ratio’s are much more practical than for an 
earth-twin, and it may not be necessary to block the star. 

Although the transit signal is of short duration (minutes), 
the combined signal can be detected continuously (for years). 
Consider the Fourier series of a rapid transit; 


r!S n 


1 - 


cos^ 6 


CO 

n=l 


Sin nojT 

-cos nojt 

n(jjT 


(A.l) 

The short duration of the eclipse gives rise to peak amplitudes of 
a fraction jrS of Iq (of the order 10“^/o). The Fourier coeffi¬ 
cients of the transit signal are 


8r sin nwT I r^cos^O 
" nS na)T \ ’ 


(A.l) 


which are of the order r/S compared to the coefficients in the 
phase signal, which are of order 1. In other words, in the Fourier 
spectrum, the reflected signal is about a factor S jr weaker (10^^) 
than the transit signal. To observe the phase part in a combined 
phase-plus-transit signal, the noise from the transit must be less 
than the phase peak amplitudes. Our Sun has noise levels of 10“^ 
at (transit-frequencies of) 10~ mz. and a noise level 10“^ for (or¬ 
bital frequencies) 10 ®Hz, see Aigrain et al. |2004 1 , so the large 
transit component does indeed pose a problem since the abso¬ 
lute noise will be large. Of course, when a planet is periodically 
blocking its star, the star is also blocking the planet every other 
half period. The corresponding dip in the signal is again a factor 
S jr smaller. 

The effect of discontinuous but periodic observation gaps 
on the Fourier peaks is different from that of a transit. For ob¬ 
servation gaps, the continuous signal is multiplied with a peri¬ 
odic function, which is zero for the intervals sin|wf| < sinwr, 
whereas for a transit, such a function is added instead. Periodic 
interruptions at the inferior and major conjunction are inevitable 
for an edge-on observation, when the planet passes behind an oc- 
culter or within the telescope’s inner working angle. The block¬ 
ing durations It must obviously be short (compared to the pe¬ 
riod). In the short-gaps limit, the multiplicative function be¬ 
comes 1 - 2 wt 2 e^‘"‘^'/7r. The even and odd peaks, / 2 „ and 
/ 2 „+i are shifted by (+/(0) - f(nla)))a>Tln, in the approxima¬ 
tion naiT « 1 . 
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Appendix B: Kepler orbits 

For planets with noncircular, but eccentric orbits, the distance to 
the host star is a time-varying periodic function. As a result, the 
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^ OJt 


so that 

- Aojt - -d/l. 

r(t) a 

When the complex harmonic e'""' is expressed as a periodic 
function of A, one obtains the generating function for the Bessel 
functions (of the first kind) 

OO 

^ gm^g-iH£sini ^ ^ 

k=-oo 

Combining the last two identities and integrating over a full pe¬ 
riod gives us the Fourier coefficients for l/r(f), which are 

TT 

^ re‘""'-^dtuf =-7„(ne). 

In J r{t) a 

-TT 

Hence, the reciprocal distance has Fourier series 


Fig.A.l. Phase-light curves for top view on an eccentric plan¬ 
etary orbit. Values for the eccentricity are e = 0 (gray), e = . 1 
(blue), .2 (black), and .3 (purple, highest peak intensity). The 
visuals and the inset are for e - .2. 

light intensity received by the planet is periodically modulated. 
We will calculate here the effect of (small) orbital eccentricity 
on phase light curves. The position vector is no longer given by 
the circle Q but with an ellipse with a semimajor axis a along 
the x-direction and eccentricity e, given by 

r(f) = x{t)i H- y{t)j - (i cos v(f) -H j sin v(f))r(f). 

The Cartesian coordinates (x,y,z) - (x,y, 0) are 

x(f) = a cos A(t) - ae, y{t) - Vl - e^a sin /1(f), 

and the polar coordinates (r, p, v) = (r, |, v) are given by 

// // cos A(t)-e 

r{t) - a- aecos /1(f), cos v(f) = -. 

1 - ecos /1(f) 


1 


1 

- y 7„(ne)e'”"' 
a ^ 

n=-cc 


1 

a 


2 

a 


CO 

^ J„(ne) cos nojt. 

n=l 


We used - Jn{ne) to get the last identity. It is in¬ 

teresting that the Fourier coefhcients of r(f) and of cos 0(f) are 
Bessel functions, although the functions r(f) and cos 0(f) them¬ 
selves have no explicit expression in f. The Fourier coefficients 
of the all powers r^(f), with V > -1 of the radial coordinate r(f) 
are finite combinations of Bessel functions too. 

The luminosity of a planet is, however, inversely propor¬ 
tional to the square of the distance r(f). We may use the con¬ 
volution the series for l/r(f) with itself to obtain 


1 

rHt) 


- CO CO 

— ^ ^ Jkike)J„^kine - ke). 

rt=-oo k=-oo 


Alternatively we may write 


1 ^^f ^ 1 i _ 

r^{t) 1-ecos/l 


-dA 


and use the series 


These orbital coordinates are expressed in the eccentric anomaly 
A (usually denoted E, see Goldstein|1964| l. This parameter A is 
(implicitly) defined as a function of the orbital phase wf, from 


tot - A - esin/l. 


An expansion in e for the coordinates in the time variable can be 
obtained with 

/1(f) - cot + e sin ^wf -H e sin (cot -!-■■■)). 


1 

1 - e cos A 


VT^ 


f b 


z 




The parameter e is a new measure of the orbital eccentricity (also 
in the interval [0,1]) defined with e = tan (2 arcsine). We may 
obtain a second way to express the Fourier series 


1 

rHt) 


n=-oo Vi € k=-oc 


In the general case, the eccentricity of the orbit leads to a 
modulation factor jr^it) in the light curve. The simplest case 
arises for a homogeneous planet with inclination 0 = 0. (Taking 
Mi = 1, e = 0, the half-illuminated disk gives f - neglecting 
the possible 45 degrees glint effect.) The signal is 


We can now pick out the Fourier coefficients of the light curve. 
They are, for n > 1, 


/« = 


16 

3n Vl - 


k~—oo 


16e" 

37r2'' 



(B.l) 


Z 6 Q 
fn COS nut = — 

„=0 ^TrrVO 


8 

37r(l - ecos/l(f))2' 


The graph appears in Fig. |A.l| Note that Mercury has e - .2. 

We calculate the Fourier coefficients fn by first considering 
the function l/r(f). We use A as integration variable and require 
for this substitution 


d(ut) - d(/l - esin/l) = (1 - ecos/l)d/l. 


The second, approximate equality contains the lowest-order con¬ 
tributions in e. The finite sum is the Taylor polynomial of degree 
n for the function exp(n). Therefore the coefficients drop to zero 
exponentially if e < 2/e, as ~ e^e" 12". 

Since Earth-like planets with a stable water-ice equilibrium 
will not likely occur for orbits that are too eccentric, we are 
mainly interested at low values of e. With the above equations 
it is possible to express the sine components as linear combina¬ 
tions in e of the cosine components for exact circular orbits. 
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